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Motivated by higher-dimensional theories that predict new 
effects, we tested the gravitational 1/r^ law at separations 
ranging down to 218 /^m using a 10-fold symmetric torsion 
pendulum and a rotating 10-fold symmetric attractor. We 
improved previous short-range constraints by up to a factor 
of 1000 and hnd no deviations from Newtonian physics. 
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It is generally assumed that non-relativistic gravity 
obeys an inverse-square law for all distances greater than 
the Planck length i?p = ^Gfijo' = 1.6 x IQ-^^ cm. How- 
ever, the gravitational interaction has only been tested 
[0-0 with good precision for separations greater than 1 
cm and, as summarized in Ref. very little is known 
about gravity at length scales below a few mm. Recently 
theorists, using several different arguments, have sug- 
gested that the unexplored short-range regime of grav- 
itation may hold profound surprises [p|-|l2|, i.e. that 
the gravitational interaction could display fundamentally 
new behavior in the mm regime. 

Many of these arguments are based on the notion, in- 
trinisic to string or M theory, of more than 3 spatial 
dimensions. To maintain consistency with a vast body of 
observations the extra dimensions must be "curled up" in 
very small regions, usually assumed to be comparable to 
i?p, or else hidden in some other way ||l^. It has recently 
been noted ||9|,[l0|] that the enormous discrepancy between 
natural mass scales of the Standard Model of particle 
physics (AfsM ~ 1 TeV) and of gravity (the Planck mass 
Afp = ^ncjG = 1.2 X 10^*^ TeV) could be ehminated 
if gravity propagates in all the space dimensions while 
the other fundamental interactions are constrained to the 
three familiar dimensions. This unification scenario re- 
quires that some of the extra dimensions have radii i?* 
that are large compared to i?p with 
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where M* is the unification scale (usually taken as Msm) 
and n is the number of large extra dimensions. The sce- 
nario with n = 1 is ruled out by astronomical data. If 
there are 2 large extra dimensions, R* must be about 
1 mm, and the gravitational inverse-square law (which 
follows from Gauss's Law in 3 spatial dimensions) will 
turn into a l/r'*-law (Gauss's Law in 5 dimensions) at 
distances much smaller than R* . 



Other theoretical considerations also suggest that new 
effects may show up at short distances; string theories 
predict scalar particles (dilatons and moduli) that gen- 
erate Yukawa interactions which could be seen in tests 
of the 1/r^ law. If supersymmetry is broken at low 
energies these scalar particles would produce mm-scale 
effects Finally, there may be some significance 

to the observation [|l2[ that the gravitational cosmolog- 
ical constant, A « 3 keV/cm^, deduced from distant 
Type lA supernovae [^5|jl^] corresponds to a length scale 
mm. These, and other, considerations sug- 
gest that the Newtonian gravitational potential should 
be replaced by a more general expression pTl] 
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The simplest scenario with 2 large extra dimensions pre- 
dicts X = R* and a = 3 or a = 4 for compactification 
on an 2-sphere or 2-torus, respectively |]l^, while dila- 
ton and moduli exchange could produce forces with a as 
large as 10^ for Yukawa ranges A ~ 0.1 mm [Tl| ]. 

This letter reports results of a test of gravita- 
tional inverse-square law at length scales well be- 
low 1 mm. Our instrument is shown in Fig. |^. 




FIG. 1. Scale drawing of the torsion pendulum and rotat- 
ing attractor. The active components are shaded. For clarity, 
we show an unrealistically large 1.5 cm vertical separation be- 
tween pendulum and attractor, and omit the BeCu membrane 
and the attractor drive mechanism. The 4 horizontal screws 
were adjusted to make the pendulum precisely level. 

An 82 cm-long 20 /im-diameter tungsten fiber, hang- 
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ing from an x-y-z-'d stage, supported a torsion pendulum 
with 10-fold rotational symmetry. The pendulum was 
placed above a 10-fold symmetric attractor that rotated 
slowly about the vertical axis of the pendulum. The ac- 
tive component of the pendulum was a 2.002 mm thick 
aluminum annulus containing ten equally-spaced 9.528 
mm diameter holes, centered on a 55.33 mm diameter 
circle. The active component of the attractor consisted 
of two coaxial copper disks. The 1.846 mm thick upper 
disk had ten 9.535 mm diameter holes centered on a 55.31 
mm circle, while the 7.825 mm thick lower disk had ten 
12.690 mm diameter holes centered on a 55.33 mm circle. 
The holes in the lower disk were rotated azimuthally by 
18 degrees compared to the upper holes; this geometry 
substantially reduced the signal from Newtonian gravity 
but had little effect on a short-range signal. 

We measured the torque on the pendulum, produced 
by the interactions between the holes in the attractor and 
the pendulum for pendulum-to-attractor vertical separa- 
tions, down to 218 ^m. Because of the 10-fold rota- 
tional symmetry of the instrument this torque varied at 
10 times the attractor rotation frequency lo. Spurious 
torques from electrical forces were minimized by plac- 
ing a tightly stretched, 11.4 cm diameter, 20 /xm thick 
beryllium-copper membrane between the attractor and 
pendulum, 48 /xm above the upper surface of the attrac- 
tor. The pendulum, mirrors, membrane and upper sur- 
face of the attractor were coated with gold, and the pen- 
dulum and membrane were surrounded by a gold-plated 
copper enclosure. Crucial dimensions of the pendulum 
and attractor were good to ±2.5 /xm, and the membrane, 
which had a mirror-like finish, was flat to within ±15 /xm. 

The torque on the pendulum was inferred from the 
pendulum twist angle, 9, which we measured by reflect- 
ing an autocoUimator light beam twice from either of 2 
flat mirrors mounted on the pendulum. The attractor 
rotation period, ta, was set to ITtq, where tq = 2n/uJo ~ 
399.6 s was the period of pendulum's free torsional os- 
cillations. The attractor angle, ((> = cut, and signals from 
the autocoUimator, the ion-pump pressure gauge, a 2- 
axis electronic level on the instrument, and 8 tempera- 
ture sensors were digitized at regular intervals, typically 
Tu = To/36, and stored in a small computer. Data, 
acquired continuously, were broken into "cuts" each of 
which had exactly 2 oscillations of the lOw signal. For 
each "cut" , the filtered ||l^] pendulum twist 9 as a, func- 
tion of (f) was fitted with 
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= 5Z sinn0 ± c„ cos ncj)] + ^ dm-Pm , (3) 

n m— 

where the first sum ran over n—l, 3, 10, 20 and 30. The 
floating n=10, 20 and 30 terms yielded the fundamen- 
tal and leading harmonics of the torque signal; the small 
n = 1 and n = 3 terms were fixed by fitting data cover- 
ing one complete revolution of the attractor. The n=3 
coefficients accounted for the 3-fold symmetry of the at- 
tractor drive system and the n=l coefficients, typically 



consistent with zero, were used for diagnostic purposes. 
The di and c?2 polynomial coefficients allowed for the slow 
unwinding of our torsion fiber (<0.2/xr/h). 

Figure g shows a typical segment of fitted autocoUima- 
tor data. 
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FIG. 2. AutocoUimator data for one complete revolution 
of the attractor, taken at ^ = 237 ^m. The curve is a fit to 
Eq.| 

The resulting n — 10, 20 and 30 torques are 

Tn = I\(bn + ic„)(w^ - n^w^ + 2i-/nuj)\ sec(^^) , (4) 

ZUJo 

where / = 128.4 g-cm^ is the calculated rotational in- 
ertia of the pendulum, the superscript tildes on the b 
and c coefficients indicate corrections for electronic time 
constants (see Ref. [||), and the damping coefficient 
7 « 5 X 10~^ s~^ is the reciprocal of the observed decay 
time of the free torsional amplitude. We measured ojq 
for each data point and found it was constant to within 
0.3%. 

Our data sets consisted of the Tio, T20 and T30 torques 
as functions of separation. Clearly, compared to Newto- 
nian gravity, a short-range l/r^-violating force will have 
both a higher harmonic content and a stronger depen- 
dence on (. Data taken at C's ranging from 218 /xm to 
10.78 mm and at several horizontal displacements x and 
y, were compared to calculations of the signals expected 
from the potential of Eq. |^. The expected Newtonian 
(Yukawa) twists were computed to a precision better than 
0.1(0.3q!) nrad while a typical measured point had a sta- 
tistical error of 8 nrad. We found analytic solutions for 
4 of the Newtonian integrals and 2 of the Yukawa inte- 
grals and evaluated the remaining integrals numerically. 
Figure ^ shows the Tig, T20 and T30 torques together 
with the Newtonian predictions. We fitted these data, 
for a series of A values, with the physics parameter a 
and 4 constrained and 2 unconstrained parameters de- 
scribing the instrument. The constrained instrumental 
parameters (the gap e — 10 ± 10 /xm between the two 
attractor disks, the ratio of the "masses" of the holes in 
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the upper and lower disks p = 0.1329 ± 0.0001, a factor 
/ = 1.000 ± 0.003 representing the uncertainties in the 
torque scale and the "mass" of the pendulum holes, and 
= 910 ±5 /im the z-micrometer reading at C = 0) were 
determined by auxilliary measurements described below. 
The unconstrained intrumental parameters a;o and ?/o 
were the centered values of x and y micrometers. In all 
cases, the best-fit values of the constrained experimental 
parameters agreed with their nominal values, the aver- 
age and maximum deviations were 1.3(t and 1.9(7, respec- 
tively. Our data are well described by Newtonian gravity 
with a per degree of freedom of 1.03 with a probability 
of 41%. The resulting contraints on new physics of the 
form given in Eq. |^ are shown in Fig. ^; scenarios with 
a > 3 are excluded at 95% confidence for A > 190 /im. 
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FIG. 3. Tio, 720 and T2.0 torques, and the lOw torque from 
the upper attractor alone. Each point contains at least 36 
individual measurements; the Icr error bars are derived from 
the scatter of the individual measurements. The upper panel 
shows the measured torques; the solid line is the Newtonian 
prediction. The Tio sign change at z = 2.2 mm is due to 
cancellation from the bottom set of holes. The lower panel 
displays the Newtonian fit residuals; the solid curve shows the 
effect on Tio of an interaction with a = 3 and A = 250 pim. 

Accurate alignment and control of the pendulum and 
attractor geometry and precise calibrations were essential 
in this experiment. 

The attractor geometry parameters were determined 
using a measuring microscope, and the "masses" of 
the pendulum and attractor holes were computed from 
their dimensions and the measured densities of the ac- 
tive components of the pendulum and attractor. The 



torque scale was calibrated gravitationally. Two small 
aluminum spheres, inserted into diametrically oppo- 
site holes in a similar pendulum annulus, gave it a 
known 522 spherical multipole moment. A Q22 at- 
tractor moment was created by two bronze spheres, 
centered on the torsion pendulum and separated by 
27.965 cm. These spheres were rotated around the 
pendulum axis by a turntable to produce a (4.010 ± 
0.001) X 10^^ dyne-cm 2uj torque on the pendulum. 
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FIG. 4. 95% confidence upper limits on 1/r -law violating 
interactions of the form given by Eq. ^ The region excluded 
by previous work p||3|j2c| lies above the heavy lines labeled 
Irvine, Moscow and Lamoreaux, respectively. The data in 
Fig. H imply the constraint shown by the heavy line labeled 
Eot-Wash. Constraints from previous experiments and the 
theoretical predictions are adapted from Ref. [^, except for 
the dilaton prediction which is from Ref. 

The z-micrometer setting at vanishing vertical 
pendulum-to-membrane separation was inferred from the 
capacitance between the pendulum and the membrane. 
(The membrane and the upper fiber attachment, which 
were normally grounded to the vacuum can, were elec- 
trically isolated for these measurements.) We fitted the 
capacitance as a function of the z-micrometer setting, see 
Fig. H, with an expression for a plane electrode parallel to 
an infinite plate [|l9| . (The function, which contained a 
leading-order expression for the fringing fields, was inte- 
grated over the pendulum's sinusoidal "bounce" motion.) 
The fixed vertical separation between the membrane and 
attractor (48 ±3 /xm) and their parallelism (misalignment 

< 0.2 mrad) were measured using depth micrometers and 
a special jig that could be substituted for the membrane. 
A capacitive measurement gave a consistent value for this 
separation. The peak-to-peak vertical runout and "wob- 
ble" as the attractor rotated were ^1.3 /im and ^ 5 /im, 
respectively. 

The pendulum plane was adjusted to be horizontal to 

< 0.2 mrad by temporarily replacing the normal mem- 
brane with 2 isolated semicircular copper disks that, to- 
gether with the pendulum, formed a differential capac- 
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itor. We then monitored the the lu differential capac- 
itance signal as d was rotated at a steady rate and ad- 
justed the pendulum's small trim screws (shown in Fig. |^) 
to minimize this signal. (The leveling precision was lim- 
ited by the ±74 /im horizontal runout of the fiber as 
■d was rotated.) The semicircular disks were removed 
and the reinstalled membrane was made parallel to the 
pendulum (to within < 0.1 mrad) by tilting the entire 
instrument and finding the angle that gave the smallest 
capacitance for a fixed value of z. Then the x and y mi- 
crometers were adjusted to recenter the pendulum on the 
attractor by maximizing the lOcj, 20a;, and 30a; signals. 
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FIG. 5. Capacitance measurements of the separation. The 
horizontal axis is the reading of the z-micrometer. The 
fit indicated that the pendulum would touch the screen at 
2 = 978 ± 3 /xm. 

We checked that likely non-gravitational sources of sys- 
tematic error had negligible effects on our results. The 
attractor drive mechanism was largely nonmagnetic; the 
amplitude of magnetic field variations at the site of the 
pendulum was < 0.6 fiG at the signal frequency. The 
measured sensitivity of the pendulum to magnetic field 
variations at the signal frequency was 57 nrad/mG, cor- 
responding to a magnetic systematic error of 0.03 nrad. 
The energy dissipated in the attractor motor warmed 
the attractor slightly. A 5-fold increase in the motor's 
temperature rise, produced by attaching heat-dissipating 
resistors to the motor housing, changed 5io and 620 by 
17±11 nrad and — 4±9 nrad respectively. We found that 
an overall 98 mK change of the instrument's temperature 
changed the equilibrium twist of the fiber by 1.17 /irad. 
During normal operation, none of the sensors on the in- 
strument had temperature variations at the lOo; signal 
frequency that exceeded 100 /iK, indicating that spuri- 
ous signals from temperature variations were ^ 1 nrad. 
When the attractor rotation period was changed from 
7to to 43ro (the normal period was 17ro), 610 and 620 
were unchanged to within 5 ± 15 nrad and 11 ± 31 nrad, 
respectively. 

Our results, interpreted in the simplest unification sce- 
nario with 2 equal large extra dimensions, imply a unifi- 
cation scale given in Eq. |l| of M* > 3.5 Tev. We are now 
preparing a second-generation experiment with a differ- 



ent pendulum and attractor that should provide higher 
precision and better sensitivity at small A. 

We thank Profs. David Kaplan and Ann Nelson for 
helpful discussions, and Nathan Collins, Angela Kopp 
and Deb Spain for assistance with the experiment. This 
work was supported primarily by the NSF (Grant PHY- 
9970987) and secondarily by the DOE. 



[1 
[2: 

[3] 

[4] 

[6; 

[7 

[8 

[9: 
[10: 

[11 

[12 
[13 

[14 
[15 
[16 
[17 

[18 
[19 

[20 



D.E. Krause and E. Fischbach, |arXiv:hep-ph 9912276 . 
J. K. Hoskins, R. D. Newman, R. Spero and J. Schultz, 
Phys Rev. D 32, 3084 (1985) 

V. P. Mitrofanov and O. I. Ponomareva, Sov. Phys. JETP 
67 1963 (1988) 

G.L. Smith, CD. Hoyle, J.H. Gundlach, E.G. Adel- 
berger, B.R. Heckel and H.E. Swanson, Phys. Rev. D 
61 022001-1 (2000). 

Y. Su, B. R. Heckel, E. G. Adelberger, J. H. Gundlach, 
M. Harris, G. L. Smith, and H. E. Swanson, Phys. Rev. 
D 50, 3614 (1994). 

S. Baefiler, B. R. Heckel, E. G. Adelberger, J. H. Gund- 
lach, U. Schmidt, and H. E. Swanson, Phys. Rev. Lett. 
83 3585 (1999). 

C.M. Will in The Confrontation between General Rel- 
ativity and Experiment, eds. M.C. Bento et al., World 
Scientific, Singapore, p. 183 (1993). 

G. C. Long, H. Chang and J. Price, Nucl. Phys. B 529, 
23 (1999). 

N. Arkani-Hamed, S. Dimopoulos and G. Dvali, Phys. 
Lett. B429, 263 (1998). 

N. Arkani-Hamed, S. Dimopoulos, G. Dvali and N. 
Kaloper, Phys. Rev. Lett. 84, 586 (2000). 
S. Dimopoulos and G. Guidice, Phys. Lett. B379, 105 
(1996). 



R. Sundrum, |arXiv:hep-th/9708329 . 

L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 



(1999); a,rXiv:hep-th 9906064. 



D.B. Kaplan and M.B. Wise, arXiv:hep-ph 0008116 



A.G. Riess et al, Astron. J. 116, 1009 (1998). 

S. Perlmutter et al, Astrop hys. J. 517, 565 (1999) . 

A. Kehagias and K. Sfetsos, |arXiv:hep-th/9905417 (Our 

Eq. 1^ is valid for extra dimensions of size R* as long as 

r > R'-) 

A simple digital filter suppressed free torsion oscillations 
as described in Ref. 

L.D. Landau, E.M. Lifshitz and L.P. Pitaevskii, Electro- 
dynamics of Continuous Media, 2nd erf., Pergamon Press 
(1984), pl8. 

S. D. Lamoreaux, Phys. Rev. Lett. 78, 5 (1998). 



4 



